Fractional derivatives of random walks: Time series with long-time memory 
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We review statistical properties of models generated by the application of a (positive and negative 
order) fractional derivative operator to a standard random walk and show that the resulting stochas- 
tic walks display slowly-decaying autocorrelation functions. The relation between these correlated 
walks and the well-known fractionally integrated autoregressive (FIGARCH) models, commonly 
used in econometric studies, is discussed. The application of correlated random walks to simulate 
empirical financial times series is considered and compared with the predictions from FIGARCH 
and the simpler FIARCH processes. A comparison with empirical data is performed. 
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I. INTRODUCTION 

In recent years, there has been a surging interest in the 
application of non-integer order differential operators to 
describe different types of temporal and spatial anoma- 
lies displayed by complex systems fEHllIllllSISHflU- 
For example, a fractional diffusion equation suitably de- 
scribes the asymptotic behavior of independent random 
walkers on random fractal structures and also corre- 
sponding fractional Fokker-Planck equations have been 
used to model anomalous behavior d • Several recent re- 
views can be found in literature covering a vast field of 
applications (see e.g. [1, @). A treatment from a math- 
ematical point of view can be found in [3, 4]. Closely 
related to the present work is the subject of fractional 
derivative operators as discussed within the framework of 
random walks @,Q and anomalous transport phenomena 

!. 

Random walks play also an essential role in model- 
ing the time evolution of stock prices. The literature in 
the field is huge and we just mention few paper s related 
with the present work [3, El El El, El El M, El El • 
The different features studied include stock-index cross- 
correlations EoKprobability distribution functions of log- 
returns [H EMI, El El, USi El El , leverage effect [lj| , 
stock-stock cross-correlations and market internal struc- 
ture G3, tick-by-tick behavior [l~7j . non-stationarity is- 
sues [l8j . 

A succesful and widely used model to describe price 
variations is based on an autoregressive process with con- 
ditional heteroskedasticity (ARCH) due to Engle (see e.g. 
[ll), of which simple variants have been recently sug- 
gested [H, HJ . Further variants of ARCH models have 
been recently discussed with regards to volatility [H, HI . 
ARCH type models have been generalized to incorporate 
a long-time memory in the surrogate time series, in an 
attempt to mimick the strong autocorrelations observed 
in volatility and absolute returns [H HE US HI HI, HI • 

The issue of long-time autocorrelations, or long time 
memory, in time series is intimately related with the con- 
cept of Hurst exponent, introduced by Hurst many years 
ago to describe the persistence observed in the behavior 
of Nile floods [30(. Motivated by these ideas, Mandelbrot 



introduced few years later a long-time memory model 
known as fractional Brownian motion (FBM) , being sta- 
tionary on all time scales, as an attempt to model such 
anomalous phenomena [3l| . Rigorous mathematical as- 
pects of FBM can be found in textbooks [32| , the latter 
including also a discussion of Levy processes with long- 
time memory. A recent work considers other aspects of 
the FBM model [H. 

The problem of accurately determining the Hurst ex- 
ponent of a time series has been the subject of intense 
activity, started with the description of persistence in 
DNA sequences in which the detrended fluctuation anal- 
ysis (DFA) has been introduced Further develop- 
ments were achieved based on Haar wavelets (35j and its 
generalization to higher orders [36| . In the present work, 
we will make use of Haar wavelet techniques to analyse 
the surrogate time series and determining the associated 
Hurst exponents. 

In this work, we apply a fractional derivation (and 
integration) operator to an uncorrelated random walk, 
obtained within a simple ARCH prescription, to gener- 
ate a surrogate financial time series with uncorrelated 
variations on all time lags, but having a slowly decay- 
ing auto-correlation for the absolute variations, repre- 
senting the absolute returns. The new process is denoted 
as fractional random walk ARCH (FRWARCH) and we 
show that it has finite second moments, in contrast to FI- 
GARCH processes characterized by an infinite variance. 
We will conclude that FRWARCH can become a useful 
tool in econometrics applications. 

The paper is organized as follows. In Sect, m we 
briefly review the fractional derivative operator and its 
main properties. In Sect. IIII1 we discuss the effects of 
applying a fractional derivative (and fractional integra- 
tion) operator to a standard random walk. The scal- 
ing properties of the resulting fractional random walk 
are discussed. Sect. |TV] is devoted to FRWARCH, for 
which the associated probability distribution functions 
and auto-correlations are studied. In Sect. El we review 
the widely used FIGARCH process and discuss several of 
its properties in comparison with those of FRWARCH. In 
Sect. I VII we briefly consider, for completeness, the sim- 
pler fractionally integrated ARCH models currently used 
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in literature. In Sect. IVIII we present a brief comparison 
of FRWARCH and FIGARCH with empirical data to see 
how they actually perform in more realistic situations. 
Reference is made to a simpler generalized ARCH model 
without long-range memory, but having a slow decaying 
autocorrelation function, to better understand the role of 
memory in the models presented here. Finally, Sect. I Villi 
summarizes the main conclusions of the work. 



II. FRACTIONAL DERIVATIVE OPERATOR 

Let us consider a time-dependent function y(t) 
recorded at times t — ir, where < i < N and r is 
the time resolution. For simplicity, we will indicate the 
associated time series as y^. 

In the following, we study the fractional derivative 
(finite difference) operator of fractional order a, Di a \ 
known in literature as the Griinwald-Letnikov scheme 
(see e.g. @, Q), which is defined as 



i=0 



where 



r(a + l) 



r(» + i)r(a + 1 - ») 



(i) 



(2) 



is the binomial coefficient and T(x) is the Gamma func- 
tion d3. 

For the numerical implementation of Eq. (jTJ) , it is more 
convenient to work out an equivalent expression for the 
binomial coefficients in Eq. @. This will also allow us 
to make contact with long-range memory models known 
from the financial literature. Using the properties of the 
r(x) function [37[, one can show that 



a\ a(a — l)(a — 2). ..(a — i + 1) 



(3) 



where the numerator is a polynomial of ith degree in a. 
Using Eq. ©, Eq. flU becomes 



r 'Vn 



= — I y« ~ aVn-i + ^(a - l)y n -2 

+ (-l)"~(a - 1)... (a -n + !)»,]. (4) 



In what follows, we will consider values of a in the 
range —1 < a < 1. Indeed, it is easy to see from Eq. (|4|) 
that the case a = 1 corresponds to the (Euler) first order 
derivative of y n , i.e., 

D^Vn = -(Vn -Vn-l), ( 5 ) 
T 

while the value a = — 1 yields the integral ofy n , 
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FIG. 1: (color online) The absolute value of the coefficients 
IC^'l vs index i, for different values of the fractional order 
a = —0.2 (squares), -0.1 (circles), 0.1 (diamonds) and 0.2 
(triangles). The symbols represent the numerically evaluated 
coefficients using Eq. (0 and the straight lines the asymptotic 
result, Eq. ©. The results for a — 0.2 have been multiplied 
(down shifted) by a factor 1/5 for clarity. 



It is also instructive to derive the asymptotic behavior of 
the coefficients in the expansion of the fractional opera- 
tor. To do this, we write Eq. © in the form 



a\ = i (0-a)(l-a)(2-a)...(i-l-a) 
i J i\ 



Now, using the relation n^oC? + z ) = r (* + z ) l T ( z ) S3 
in Eq. ©, with z = -a, the coefficient c\ a) = (~f) l (") 
can be written as 



c (a>) = T(i + z) 



r(z) *! 

T(i — a) c 
T(-a)T(i + l) ~ ~r(l-a) T(i + 1) 



T(i - a) 



(8) 



where the known relations T(l + z) = zT(z) and i\ — 
T(i + 1) have been used. Exploiting the large argument 
expansion T(az + b) ~ exp(— az){az) az+b ~ 1 / 2 , one can 

obtain the asymptotic behavior of c\ a ^ for large i as 



C 



(a) 



,-(H 



r(i 



i > I, 



(9) 



which is positive for a < and negative for a > 0. Re- 
sults for |C| Q) | are plotted in Fig. Q] in double logarithmic 
scale for few representative values of a. To be noted is 
the rather quick approach of the actual value of C[°^ to 
its asymptotic form Eq. ([9]). This feature can be used 
in numerical calculations where c[ a ^ can be replaced by 
Eq. ([9]) when the difference between the two is sufficiently 
" small. 

£>[ -1) 2/n = T(y„+y„_i+y„_2H hyo) = X^ r - ^ lt is natural to refer to as a fractional derivative 

i=o operator for values < a < 1 and as a fractional integral 
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FIG. 2: (color online) The fractional function yn' = D±*'y n 
vs index n, in the case y„ = 1 (i.e. 7 = 0), r = 1, for different 
values of the fractional order: a = 0.4 (squares), 0.2 (circles), 
-0.2 (diamonds) and -0.4 (triangles). The symbols represent 
the numerical results using Eq. , and the straight lines the 
asymptotic form, Eq. (Ilip . The dashed line corresponds to 
the original series y n . 



operator when — 1 < a < 0. This can be seen in the 
case that y n ~ n 7 . In fact, the action of the fractional 
operator, which will be denoted as 



Vn = V 



(a) 



1 " 

-yc 



(«)„ 



(10) 



i=0 



can be estimated for large n using the asymptotic be- 
havior for c[ a \ Eq. ©. Roughly, we can write y« ' — 
f™dx x- {1+a ^ (n - a;) 7 . When 7 > and n > 1, the 
largest contribution to the integral comes from the small- 
est values of x, yielding y^ ~ n 7 J™ c?x 2; -( 1 + Q; ) ~ n 7_Q . 
Thus, 



(ii) 



We show in Fig. [5] numerical examples of the fractional 
operation D^y n (i.e. r = 1) in the case y n — 1 (7 = 0), 
as an illustration of the result Eq. (fTTj) . Note that for a 
constant function y n , the slope of the curves is just minus 
the order of the fractional operator. To be noted also is 
the rather quick approach of the fractional operator re- 
sult Eq. (flO|) to its asymptotic behavior Eq. (fTTj) , already 
for small values of n. 
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FIG. 3: (color online) Fractional random walks: Plots of y„ 
versus time step n, for selected values of a as indicated on the 
plot. The standard RW corresponding to a — 0, y^ = y n , is 
also shown for comparison. The walks consist of N = 32000 
steps each. 



and unit standard deviation, i.e. (77) = and (?y 2 ) = 1, 
yielding a v — 1. Thus, 



Vn 



i=l 



1 < n < N, 



(12) 



where we take yo = for simplicity. We are interested 
in studying the behavior of the associated fractional ran- 
dom walks (FRW's), obtained by applying the fractional 

operator to the random walk 'profile' y n , Eq. (|12|) . 

D {a) v - ?; (q) 

Examples of fractional walks y„ are illustrated in 
Fig. [3] As is apparent from the figure, the amplitude 
of the fractional walks increases for negative a with re- 
spect to the uncorrelated case y„ , while for positive a 
the amplitude of the walks gets smaller. The question 
is whether the fractional walks are also statistically dif- 
ferent from the uncorrelated case in the sense that long- 
range autocorrelations are present. Before discussing this 
quest, let us consider next the distribution function of the 
increments of a FRW. 



A. Probability distribution functions 



III. FRACTIONAL RANDOM WALKS 

Let us consider the case in which y n represents a stan- 
dard random walk (RW), and assume the time resolution 
t = 1 for convenience. Specifically, y n is constructed as 
a sum of independent, equally distributed random num- 
bers, r). Without loss of generality, let us assume the lat- 
ter are drawn from a Gaussian distribution of zero mean 



The increments of a fractional random walk are de- 
noted as 



(a) 



Vn a) 



(a) 
Vn-V 



(13) 



For a standard random walk, as in Eq. (|12p . one simply 
has Ay n = rj n , which is a local function of the time step 
n, i.e. independent of the previous values of rji, i < n. In 
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contrast, for any fractional order a 0), we find 



n-l 



l-i 



i=0 



8=0 

n-l 



(14) 



where we have taken yo = 0. Thus, for fractional a, the 
increments of the fractional random walk, Ay n a ^ , require 
the knowledge of the whole history of the walk, {rji}, with 
1 < i < n. This result suggests that FRW's may possess 
long-range autocorrelations, as we will indeed see below. 

The probability distribution function (PDF) of the in- 
crements, denoted as G(g), is a function of the scaled 

variable g = g n = {Ay n a) - (Ay n a ^))/a a , where a a is 

the standard deviation of Ay n a \ If r\ is normally dis- 
tributed, G(rj) = G {r,) = (Wp 1 / 2 exp(-7y 2 /2), so is 
also G(g). This does not hold in general for distributions 

G(ri) ^ G (77) because the coefficients Cj in Eq. (fl"4|) . 
weighting the r/n-i's differently, can yield effectively non- 
indentically distributed random numbers and the central 
limit theorem does not hold. To see this in a simple case, 
consider that = exp(— ija). Now, if a — > then 

G(g) = G(rj). Also for finite a <C 1, the distribution 
G(g) will not be Gaussian. In the present case, we illus- 
trate this behavior by assuming that G{rf) = exp(— \rj\). 
The corresponding PDF's for g are shown in Fig. Q] in 
the case a = —0.2 for series of different length N, 
where < n < N. As one can see from the figu- 
re, G(g) remains exponentially distributed attaining the 
form G(g) ~ A exp(-\g\/A), with A ~ 0.7. 



B. Auto-correlations 

Since for a standard RW the profile behaves, in a 
statistical sense, as y n ~ n 1 / 2 , actually meaning that 
({Un — yn+m) 2 ) ~ \ m \i the fractional random walk obeys 



(15) 



corresponding to ((j/n — y^+ m ) 2 ) ~ l^l 1 2a - Identify- 
ing in Eq. (fT5|) the power of n with the Hurst exponent 
if [H we find, 
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FIG. 4: (color online) PDF, G(g), for fractional random walks 
variations Ay^ as a function of the scaled variable g = 
(Ay ( n a) - (Ay ( n a) ))/a a in the case a = -0.2 (see Fig.©, for 
series of different length N (0 < n < N) indicated in the in- 
set. The continuous line is the normal distribution, while the 
dashed line is the exponential form f(g) — Ai exp(—\g\/A2), 
with A\ = 0.72 and Ai = 0.7. 

the auto-correlation function of the increments Ay^\ 
C a {m), which is given by C a (m) = {g n g n+m ). If 
long-range memory is indeed present, the autocorrela- 
tion function is expected to obey the scaling behavior 
C a {m) - |m|- 2 ( x - H ) (see e.g. [Mill), yielding 

C a (m) ~ -sign(a) |m|- (1+2Q) , \m\ > 1. (17) 

To detect the presence of long-range memory for a frac- 
tional random walk y^ , we apply the method known in 
literature as the fluctuation analysis (FA) [36( based on 
Haar wavelets (HW) [35j , consisting in studying the scal- 
ing behavior of y^ on the time scale t. To do this, the 
total number of points in the series, TV, is divided into 
consecutive non-overlaping segments of length £ > 1, cor- 
responding to the time scale t = It. Inside each segment 
s, 1 < s < N/£, the average of y^\ denoted as B s (£), is 
evaluated according to 



(18) 



i=i 



The FAHW approach consists in studying the fluctua- 
tions of the profile on the time scale £ = t/r, defined as 



H 



a, 



and 



H. 



(16) 



Here, we are interested in cases where < H < 1, yield- 
ing for a the range of variation — 1/2 < a < 1/2. The 
case H = 1/2 corresponds to standard (uncorrelated) 
behavior, values H > 1/2 indicate persistence or long- 
time autocorrelations, while values H < 1/2 yield anti- 
persistence or negative autocorrelations. 

The Hurst exponent also determines the behavior of 



F(t) = {[B s {£)-B s _ 1 {£)Y) 



2\ 1/2 



(19) 



corresponding to the first-order Haar wavelet, and the 
average is performed over all consecutive boxes s and 
s — 1. To be noted is that Eq. (fl3|) can be generalized to 
higher order wavelets (3o| . 

The dependence of F(£) on £ is expected to obey a 
scaling behavior of the form 

F{£) - £ H , (20) 
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FIG. 5: (color online) Fluctuation analysis of fractional ran- 
dom walks: Plotted is the fluctuation function F(£)/F(l) vs 
time scale i for different values of fractional order: a — 0.2 
(triangles), 0.1 (diamonds), -0.1 (circles) and -0.2 (squares). 
The straight lines have slopes H = 1/2 — a given by Eq. (I16|l . 
The series considered consist of N = 10 time steps each. 

from which one can obtain the Hurst exponent H. Re- 
sults for the fractional orders considered in Fig. [3] are 
reported in Fig. [5] The results support the expectations 
of the existence of long-range memory for FRW's. 

IV. FRACTIONAL RANDOM WALKS FOR 
FINANCIAL TIME SERIES 

It is well known that daily log-returns, i.e. the vari- 
ation of the logarithm of asset price P n at day n, de- 
noted as AS n = lnP„ — lnP n _i, is not correlated with 
its variations AS n '=£ n and the corresponding PDF dis- 
plays power-law tails, G(g) ~ \g\~~^ with j3 w 4 (see 
e -g- In addition, absolute log-returns, that is 

I ASn I , appear to be long-range auto-correlated (see e.g. 
[PA |38|). We aim at modeling these features using the 
FRW described above for negative values of a. 

In order to obtain a PDF for log-returns displaying 
power-law tails, we resort to a simple auto-regressive con- 
ditional heteroskedastic model (ARCH) [Ijj]. We define 
the log-returns according to 

AS n = CTnVn, Tl > 1, (21) 

where rj n are uncorrelated random numbers drawn from a 
normal distribution with zero mean and unit variance, i.e. 
(r)l) = 1 and (AS n AS m ) = {(AS) 2 )S n , m - The standard 
deviation a n changes in time according to the ARCH 
prescription, 

a 2 n = a + b(AX n a \) 2 , (22) 
where AX n a \ is proportional to the fractional integral of 

ASn-i, Eq. dinD, 

= , -, M a) AS n - U a < 0, (23) 
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FIG. 6: (color online) The normalization coefficients A n (a) 
versus n for different values of a (from top to bottom): a = 
-0.4, -0.37, -0.35, -0.30, -0.20 and -0.10. Notice the slow 
convergence of the coefficient for values a — > —0.5. 



and the factor A n -i(a) — y J2"=o i^i"^ ) 2 ensure s the 

constancy of the 2nd moment ((AX^"\) 2 ). Indeed, using 
Eq. (fTUl) the latter becomes 

^((AS) 2 )=al (24) 

where we find a\ = a/ (1-6) according to Eq. ^ . Plots 
of A n (a) are shown in Fig. [6] The present combination 
of a FRW with ARCH will be denoted for brevity FR- 
WARCH. 

In the following, we discuss numerical results to illus- 
trate the implementation of FRWARCH. We consider the 
fractional order a = —0.4, together with the ARCH pa- 
rameters a = 0.5 and b = 0.7. To improve the accuracy 
of the numerically obtained PDF of AS n , we have consi- 
dered walks of up to TV = 10 5 time steps and averaged 
over 100 configurations. The results for the PDF, G(g), 
are plotted in Fig. [7] as a function of the scaled vari- 
able g. A simple fit has been conducted to the numeri- 
cal results suggesting a broad distribution function with 
power-law tails, G(g) ~ \g\~^, for \g\ 3> 1, with f3 ~ 3.7, 
the latter consistent with the value expected analytically 

The fluctuation analysis of FRWARCH time series is 
reported in Fig. where we plot the function F{t) as a 
function of time scale I for AS n and its absolute value 
|A5„|, together with the behavior of the fractional pro- 
cess AX n a \ Eq. (|2"3"]) . As expected, log-returns AS n dis- 
play uncorrelated behavior (Has — 0.5), while the corre- 
sponding absolute returns show persistence on long time 
scales with an effective exponent £?|as| — 0-8- 

We also include the results for AX n a ^ indicating that 
for the fractional process the effective Hurst exponent, 
Hax — 0.87, is close to the expected value H = l/2—a = 
0.9. We conclude that the long-time memory induced in 
absolute returns is a bit weaker than the one for the frac- 
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FIG. 7: (color online) PDF, G(g), for FRWARCH varia- 
tions AS n as a function of the scaled variable g — (AS n — 
(AS n ))/a A , in the case a = 0.5, b = 0.7 and a = —0.4. The 
surrogate series considered consist of N = 10 5 time steps each 
and averages over 100 configurations were performed. The 
continuous line is a fit with the form: / = /o/(l + Ig/go] 13 ), 
with /o = 0.51, go = 0.86 and (3 — 3.7. Here, we find 
a A ~ 1.30, in agreement with Eq. (124J1 . 



O AX W (H=0.87) 
□ IAS J (H=0.80) 
O AS° (H=0.50) 




FIG. 8: (color online) Fluctuation analysis for FRWARCH 
time series in the case a = —0.4 (see Fig. [7J). Plotted is 
the fluctuation function F(£)/F(l) vs time scale £ for: Log- 
returns (circles), absolute log-returns (squares) and fractional 
random process, Eq. (|23p . (diamonds). The straight lines have 
slopes H as indicated in the inset. The fluctuation analysis 
(FAHW) was done on a single time series of length TV = 10 5 . 



tional process AXn^ ■ However, the value of -H)as| — 0.8 
is still consistent with the time decay of the autocorre- 
lation function observed in empirical data for which the 
corresponding Hurst exponent typically varies within the 
range 0.7 < If | AS , < 0.9 0, H| . 



V. FIGARCH REVISITED 

We briefly review the main properties of FIGARCH 
models [27j relevant for our present discussion. Let us 
start with the GARCH model (25|, which is considered 
here in its simplest form. It is defined analogously to 



Eqs. (J2TJ) and ([22]), where now 

a 2 n = a + b (AS n -i) 2 + c er 2 ^, 



(25) 



and a, 6, c are positive constants. The mean variance be- 
comes 

o a 



(J, 



G 



1 - (b + c) ' 



(26) 



which is finite provided that b + c < 1. In order to dis- 
cuss the FIGARCH model, one resorts to the fractional 
differencing operator, defined as 

(i-Lf = i-eL-±e(i-9)L 2 



3! 



(1 -9)(2-0)L 3 



1-][>(0)L\ 



(27) 



i=i 



where < 8 < 1, C,-(0) > and L is the lag operator, de- 
fined according to I'a;^ = x 2 l _ i . Analogously to Eqs. ([3]), 
© and the coefficients Q(6>) above behave for large 
i according to 



Ci(0) 



r(i- 



,-(!+«) 



i > 1. 



(28) 



Note that Q((9) = — C\ a \ for a > 0, as can be seen from 
Eq. (HI). Finally, taking L = 1 in Eq. (27}, one finds 
the general sum-rule, R(0) = X)i=i = 1, valid for 

< 6 < 1. 

The standard way of introducing FIGARCH is to write 
Eq. (|25|) using the lag operator L in the form 

{l-bL-cL){AS n f = a+{l-cL)[(AS n ) 2 -a 2 n l (29) 

and inserting the differencing operator (1 — L) 6 in the 
left-hand side of the above relation, yielding, 

(1 - bL - cL)(l - L) e {AS n ) 2 = a+ (1 -cL) [{AS n f - cr 2 n ). 

This is a generalization of the 9 — 1 integrated GARCH 
(or IGARCH) model [27j to the case of a fractional ex- 
ponent < 8 < 1. 

Expanding the operator (1 — L) e according to Eq. (|27p . 
we find 

al = a + b (ASn-i) 2 + c a 2 _ 1 (30) 

OO 

+ ]TC^) [(A5„_ i ) 2 -(6 + c)(AS n _ 1 _ i ) 2 ] , 



which, using the relation Ci-i(9)/Ci(8) = — 1 — 
becomes 



a +{b + 8) (AS n _!) 2 + c a 2 n _ x 

OO 

]TC^) fi(9,b,c) (AS n -i) 2 , 



(31) 



where 



fifth A (*-l)(l-ft-c)-(fl + & + c) 
fi{9, b, c) = : . (32) 



i- 1 - 
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Stability of the process is achieved when the coefficients 
fi(9, b, c), Eq. (|32[) . are positive yielding the minimal con- 
dition 



e + 2{b + c ) < 1. 



(33) 



From Eq. (|30[) . we can write the mean variance for 
F1GARCH as 



4 = 



[l-(6 + c)] [l-R(9)Y 



(34) 



which diverges for all < 9 < 1, because R{9) = 1. In 
the case < 9 < 1/2, the autocorrelation function of <r^ 
for a FIGARCH process decays as the power-law |27j 



C(r) 



-(1-29) 



T » 1. 



(35) 



(see above 



Writing the latter as C(r) ~ |t|~ 2(1 ~ 
Eq. ((ni)), we find the relation 9 = H - 1/2. 

In practical calculations, the sum in Eq. (|31j) has a 
finite number of terms, with an upper memory cut-off 
that we denote as Mq. For finite Mq, the sum R(9) < 1 
and as a result the mean variance of the process is finite. 
In order to improve the convergence of the finite M case 
to FIGARCH, one can use renormalized coefficients [27t 



Ci(0) 



so that R(9) = 
2 then reads, 



Ci(9) 



= Y,Z°iCi(9) 



(36) 



1. The final expression for 



a 



o* = a + (b + 9) {AS n -!) 2 + c <j\_ r 



(37) 



^Ci(9) fi(9,b,c) (AS n . 



if, 



i=1 



where = 6/ EZ\ C i( 9 )- 

In order to make contact with the FRW results dis- 
cussed in Sect. IIV1 we consider the simpler case c = in 
Eq. (|3"Tj) , and take a = 0.01 (in order to obtain a mean 
standard deviation similar to FRWARCH), b = 0.33 (in 
order to get a power-law PDF as in Fig. [7]) and 9 = 0.3. 
The latter is chosen to yield the same value of H = 0.8 
characterizing the autocorrelation of absolute returns ob- 
tained from FRWARCH in the case a = —0.4. Note also 
that the chosen value of b obeys the boundary condition 

Eq. U2SD. 

Results for the PDF of AS n , obtained in the case 
N = Mq = 10 5 , are shown in scaled form in Fig. [9] Also 
here, G(g) ~ Iffl - ' 3 , for \g\ 3> 1, with /3 ~ 3.9, similar to 
FRWARCH in Fig. Finally, results of the fluctuation 
analysis are displayed in Fig. I1Q[ supporting the expec- 
tation that for log-returns no correlations are present, 
while for absolute log-returns the effective Hurst expo- 
nent is consistent with H = 1/2 + 9 = 0.8. 
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FIG. 9: (color online) PDF, G(g), for FIGARCH varia- 
tions AS n as a function of the scaled variable g = (AS n — 
(AS n ))/a F , in the case a = 0.01, b = 0.33 and 6 = 0.3. 
The surrogate series considered consist of TV = 10 5 time steps 
each, with a memory cut-off of Mo = 10 time steps. (Note 
that each time series is actually 2N long, so that each of 
the N time steps considered in the average process is a func- 
tion of its previous iV steps.) Averages over 100 configu- 
rations were performed. The continuous line has the form: 
/ = /o/(l + l3/5o| /3 ), with f = 0.42, g = 1.08 and (3 = 3.9. 
Here, we find g p ~ 1.20 and 6 = 0.307. 
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FIG. 10: (color online) Fluctuation analysis for FIGARCH 
time series in the case 9 = 0.3. Plotted is the fluctuation 
function F(£)/F(l) vs time scale I for: Log- returns (circles) 
and absolute log-returns (squares). The straight lines have 
slopes H as indicated in the inset. The fluctuation analysis 
(FAHW) was done on a single time series of length N = 10 5 



with Mo = 10 steps. 



VI. FIARCH REVISITED 

There are other models currently used in literature re- 
garding long-time memory of absolute log-returns and 



variable volatility. We briefly comment on them in or- 
der to better assess the impact of FRWARCH suggested 
here. Let us consider fractionally integrated ARCH pro- 
cess [40l] . variants of which have been studied recently 
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FIG. 11: (color online) PDF, G(g), for FIARCH varia- 
tions AS n as a function of the scaled variable g — (AS n — 
(A5'n))/(7 F , in the case <to = 1-2 and 8 — 0.3 (full circles). 
The FIARCH series consisted of N = 10 5 time steps each, 
with a memory cut-off of Mo = 10 5 time steps. Averages over 
100 configurations were performed. The continuous line has 
the form: / = / /(l + \gfgof), with f = 0.39, g = 1.22 
and P = 5.7. Here, we find (vl) 1 ' 2 ~ 1.30 and 9 = 0.307. 
For comparison, results from FIGARCH (open circles, taken 
from Fig. [9} have been included. 



FIG. 12: (color online) PDF, G(g), for modified FIARCH 
variations ASn, as a function of the scaled variable g = (AS„ — 
(AS n ))/a F , in the case a = 0.67, 6 = 0.38 and 6 = 0.3 (full 
circles). The FIARCH series consisted of TV = 10 s time steps 
each, with a memory cut-off of Mo = 10 5 time steps. Averages 
over 100 configurations were performed. The dashed line has 
the form: / = / /(l + \g/gof), with /„ = 0.51, g = 0.75 
and P = 3.7. Here, we find (a*} 1/2 ~ 1.30 and = 0.307. 
For comparison, results from FRWARCH (open circles and 
continuous line, taken from Fig. [7} have been included. 



lot |41| . The model is based on the equations 
AS n = <7 n rj n , n > 1, 



CO 



Mo 

E 

i=l 



Ci(6) 



\AS„ 



\AS n \ 



(38) 
(39) 



where rj n are independent normally distributed random 
numbers = 1) and < 6 < 1/2 as for FIGARCH. 

According to the sum rule J2i=i = 1, see Eq. ([55]) . 

Eq. (|39|) yields (er„) = oo- Thus, log-returns AS n are 
uncorrelated to each other while absolute returns, |A5„|, 
display long-time memory. As one can see from Eq. (|39p . 
the only free parameters are 9 and Co [42| . in contrast 
to the three parameters for FRWARCH. This has conse- 
quences on the shape of the PDF as we can see from the 
numerical results reported in Fig. 1111 

The resulting PDF for FIARCH is indeed consistent 
with a power-law distribution at the tails, but the power- 
law exponent turns out to be large, here f3 « 6, and can 
not be controlled by tuning a model parameter. Indeed, 
the value of 9 = 0.3 is fixed by the condition that H = 
0.8. Results of fluctuation analysis for FIARCH (not 
shown here) confirm that in this case H « 0.8. Thus, it 
appears that additional degrees of freedom, in terms of 
model parameters, are required in Eq. (|39[) in order to 
obtain PDF with varying shapes. 

To do this, we suggest a slight generalization of 
Eq. (03) to the form, 



Mo 



a n = a Q J2C t (9) 



\ASn- 

(\AS n 



&|AS„-i|, (40) 



as, (er„) = cr + b (|A5 n -i|) = (To + b (a n -i)(\r)n-i\), 
yielding 



(a n ) = - 



b(\Vn\Y 



(41) 



where the mean standard deviation now can be obtained 



Note that for a normal distribution (|r?n|) = ^/2/7r w 0.8. 

Although the additional parameter b actually helps in 
getting different decaying power-law exponents (3 for the 
PDF, the model behaves similarly as FIGARCH in the 
sense that there is no an a priori way to estimate b and 
this has to be done case by case. The numerical exam- 
ples investigated (see Fig. [T2l suggest that appropriate 
set of parameters can be found but at the expense of a 
numerical search. 



VII. COMPARISON TO EMPIRICAL DATA 

In this section, we discuss results for FRWARCH and 
FIGARCH with regards to empirical data. For the lat- 
ter, we take daily data for the Dow Jones Index, Iuj{n), 
from 1st October 1928 till 16 May 2008 [H, consist- 
ing of 19994 values. As the working variable, we con- 
sider the daily changes of the logarithm of the index, 
AS n = In /dj (n) — In Idj (n— 1)- The corresponding PDF 
is shown by the full circles in Fig. [12] To be noted is that 
non-stationarity issues may play a role for the Dow Jones 
Index as discussed in [44J. Here, we disregard such cor- 
rections and assume the series as stationary. 

In the following we consider series of JV = 20000 val- 
ues. Regarding FRWARCH, we use: a = 0.5, b = 0.7 
and a — —0.4. FIGARCH results are obtained for the 
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FIG. 13: (color online) PDF, G(g), for variations AS n as a 
function of the scaled variable g = (AS n — (AS n })/cr for: 
Dow Jones Index (full circles), FRWARCH (open squares), 
FIGARCH (filled diamonds), and GARCH (open triangles 
down). The model parameters are indicated in the text. 
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time scale 

FIG. 14: (color online) Fluctuation analysis versus time scale 
[days] of absolute returns, for: Dow Jones Index (full circles), 
FRWARCH (open squares), FIGARCH (filled diamonds), and 
GARCH (open triangles down). The straight lines have slopes 
H = 0.9 and H = 0.5, and are displayed for illustration. 



case: a = 0.015, b = 0.2 and 6 = 0.4 (here also c = 0). In 
addition to these long-memory models, we consider for 
illustration a short-memory model such as the GARCH 
version discussed in Eq. (f2"5"|) , for the model parameters: 
a = 0.2, b = 0.09 and c = 0.9. Alternatively, we will 
also discuss the set b = 0.9 and c = 0.09. Results for the 
PDF's are shown in Fig. [151 To be noted is that both FR- 
WARCH and FIGARCH yield PDF's in good agreement 
with the empirical one. The results from GARCH are not 
that satisfactory, for the chosen set of paramters. If we 
take instead the alternative set (b = 0.9 and c = 0.09), 
the agreement becomes comparable to that from FR- 
WARCH and FIGARCH. The reason for the first choice 
is based on the behavior of the fluctuation function for 
absolute returns, as we see below. 

The fluctuation function of absolute log-returns for the 
Dow Jones Index is displayed by the full circles in Fig. [14] 
We find power-law behavior over about three decades 
with a Hurst exponent H ~ 0.9, indicating a strong per- 
sistence in the autocorrelation function of absolute log- 
returns. Similar behavior is displayed by lASVil for FR- 
WARCH (open squares), where the latter almost overlap 
with the empirical ones. Also FIGARCH yields results 
in good agreement with the Dow Jones values, although 
the agreement is not as good as for FRWARCH. In con- 
trast, GARCH yields good results only for time scales 
below about 100 days. On larger times, the GARCH 
fluctuation function displays uncorrelated behavior with 
exponent H ~ 1/2, indicating the existence of only short- 
range autocorrelations in the time series, as expected. 
If we use the second set of values (i.e. b = 0.9 and 
c = 0.09) for GARCH, the PDF improves considerably, 
but the fluctuation function crosses over the uncorrelated 
regime (H — 1/2) already at about 10 days, yielding a 
poor fluctuation behavior. These results suggest that a 
'long-range' memory is a necessary ingredient in a sur- 
rogate model, as the one described by FRWARCH or 



FIGARCH. 



VIII. CONCLUSIONS 

Fractional derivatives represent a conceptually simple 
scheme that allow us to build, from a standard random 
walk, stochastic processes with long-range autocorrela- 
tions. The long-memory built into the walks is controlled 
by the order of the fractional operator. Positive orders 
correspond to fractional derivatives and negative ones to 
fractional integrations. The former lead to fractional ran- 
dom walks (FRW) with negative (or anti-persistent) au- 
tocorrelations and the latter to positive (or persistent) 
autocorrelations. Long-time autocorrelations decay as a 
power-law for long time lags, the exponent of which de- 
pends on the Hurst exponent associated to the walks, 
the latter being a function of the fractional operator or- 
der. Examples have been studied to illustrate the use of 
fluctuation analysis based on Haar wavelets to determine 
the corresponding Hurst exponents. The results indi- 
cate that a constant Hurst exponent is consistent with a 
wide range of time scales, suggesting that FRW are es- 
sentially stationary for most practical purposes. A FRW 
model has been discussed for describing the behavior of 
both log-returns and their absolute values observed in 
empirical data of financial assets, which is based on a 
simple autoregressive (ARCH) scheme and denoted as 
FRWARCH. Absolute log-returns, as well as volatility, 
display strong autocorrelations, and the proposed FRW 
model seems to capture the essential features. Statistical 
properties of the present model have been compared with 
the predictions of a FIGARCH and FIARCH processes, 
to illustrate the difficulties that are found in practical 
calculations. We may conclude that FRWARCH turns 
out to be as accurate as FIGARCH regarding long-time 
memory features and it appears to be more stable than 
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the latter with regard to distribution functions of log- 
returns. We therefore suggest that FRWARCH is suitable 
for simulating empirically observed slowly-decaying ab- 
solute log-returns autocorrelations, competing with the 



presently available models in the financial literature, as 
also demonstrated by a direct confrontation with daily 
close data from the Dow Jones Index. 
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